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Abstract 

A microscopic spin model is proposed for the phenomenological Zimm-Bragg model for the helix- 
coil transition in biopolymers. This model is shown to provide the same thermophysical properties 
of the original Zimm-Bragg model and it allows a very convenient framework to compute statistical 
quantities. Physical origins of this spin model are made transparent by an exact mapping into a 
one-dimensional Ising model with an external field. However, the dependence on temperature of 
the reduced external field turns out to differ from the standard one-dimensional Ising model and 
hence it gives rise to different thermophysical properties, despite the exact mapping connecting 
them. We discuss how this point has been frequently overlooked in the recent literature. 
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Statistical descriptions of polypeptide chain conformations involve important coarse- 
graining on the level of the C a atoms. Due to planarity of the amide group, torsional 
angles of successive repeated units can be considered independent, and a pair of 0;,^ an- 
gles can be associated with each repeated unit (see Fig. Ufa)). Introduction of virtual bonds, 
connecting neighboring C a 's, strongly simplifies the description. Thus the configuration of a 
polypeptide chain can be described with the array of bond vectors {U}, i = 1...N — 1, related 
to its backbone, and a new variable 7, representing the state of the particular i— th units 
(Fig. []Jb))[l, 2|. Within the framework of the helix-coil transition theory, this variable could 
have different values in the helix and coil states, taken according to some particular prescrip- 
tion. Since Doty and co-workers experimentally showed that polypeptide chains in solution 
can be reversibly converted from the random coil to a-helix conformations [3J, a number of 
methods have been proposed to model the phenomenon [4j ranging from highly sophisticated 
computer simulations [5| to highly simplified spin-like models |6j. In the past few decades 
significant advances have been made in computational capabilities of both computer hard- 
ware and software thus enabling investigations to an unprecedented level of complexity. At 
the same time, interpretations of results from single-molecule techniques, such as stretching 
with optical tweezers g have large.y relied oa applications of dassica. spin models & In 

addition these spin approaches remain attractive for describing folding of helical proteins [91 

I I 

and influences of solvent on secondary structure formation and stability [10]. Among these, 
the Zimm-Bragg (ZB) model stands out for its success ju]. While very useful in interpreting 
experimental results, the original ZB theory is far less satisfactory from a theoretical point 
of view as it lacks a well defined microscopic description, thus preventing a clear connection 
with more sophisticated levels of theoretical description. 

It has been frequently stated in the literature, that the ZB model can be described, at 
a microscopic level, by a one-dimensional Ising model (see e.g. Ref. 12j for recent represen- 
tative examples). However, care must be exercised in making this assertion. The aim of 
this short note is directed toward formally addressing the actual equivalence of the ZB and 
one- dimensional Ising models. 

The comparative analysis will be performed by introducing a one-dimensional Potts-like 
spin model and demonstrating that it gives rise to the same thermodynamic properties as 
the original ZB model. The Potts-like model is then mapped, via an exact transformation, 
into a one-dimensional Ising-like model with an external field. It is shown that due to the 
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temperature dependence in the external field, predicted thermophysical properties are not 
equivalent to the standard one-dimensional Ising model. 

In its simplest formulation the ZB model is based on the following combinatorial rules 



111 ] . (1) Every repeated unit exists in either the hydrogen-bonded (helical) or un-bonded 
(coil) state; (2) Every un-bonded repeated unit contributes a statistical weight of unity to 
the partition function; (3) Every bonded repeated unit that follows another bonded repeated 
unit, contributes a statistical weight, s; (4) Every bonded repeated unit that follows two or 
more un-bonded repeated units, contributes a statistical weight of sa (in conjunction with 
Rule 3, this defines a); (5) Every bonded repeated unit that follows less than two un-bonded 
repeated units, contributes a statistical weight of zero. 

Note that s has the meaning of a statistical weight, and is usually interpreted in terms of 
a free energy change s = e~ l3 ^ GheUx ~ Gcoa ^ between the helix and coil states. Similar interpre- 
tation can be given for o with the additional restriction of being purely entropic in nature 



(unlike s which has both enthalpic and entropic contributions) 13] . Using combinatorial 



techniques, the ZB model allows derivation of the thermodynamic properties of the system 



from the eigenvalues of the following 2x2 matrix 11], 



M ZB = I 1 1 I (1) 



as s 



The corresponding secular equation providing the eigenvalues is, 

X 2 -X(s + l) + s(l-a) = 0, (2) 
which defines the thermodynamics of the ZB model. Note that the matrix in Eq. (PQ) is 



not symmetric. This contrasts with the results obtained 



rom the standard one-dimensional 
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As a result, eigenvalues, and 



Ising model with an external field, which is symmetric 
hence the thermodynamics are different in the two cases. We shall return to this point 
later on where comparisons between the ZB model and one-dimensional Ising model with a 
particular external field will be discussed. 

Our microscopic formulation of the ZB model is founded on a Potts-like formulation akin 



to the more general model, previously presented within a slightly different context [15 ]. 
Assume that spin ji describing the state of the z-th repeated unit can take one of Q(> 2) 
values; 7» = 1 corresponding to values of the torsional angles (pi and ipi from the helical region 
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of the Ramachandran map, while the other Q — 1 identical values correspond to torsional 
angles from the coil region. As we shall see, the condition Q > 2 plays a fundamental role 
as it accounts for the large degeneracy (and hence entropy) of the coil state. The magnitude 
of Q can be identified with the ratio of the allowed region area versus helical region area 
on a Ramachandran map. Although the formation of one hydrogen bond fixes the values 
of three couples of torsional angles (see, e.g. [lsj] ) , we simplify the model and consider only 
a nearest-neighbor construction. The energy of interaction is assumed to be different from 
zero when both 7, and 7j + i are equal to 1. The corresponding spin Hamiltonian is, 

N 

-/3 J ff = /3J^(J( 7i ,l)«J( T<+lj l), (3) 



i=l 



where 5's are Kronecker symbols and = 1/kbT is the inverse thermal energy. Here and 
in the following, we consider open boundary conditions with the iV — > 00 limit taken at the 
end. We further remark that the Potts-like model defined in Eq. flS]) differs from the classical 
Potts model proposed by Goldstein 16[ for the helix-coil transition (a detailed comparison 
between the two models is interesting and will be addressed in future work). 

The partition function Z can be obtained via standard transfer matrix techniques 

z = yj e -^(,,}, = yjn(M) w , (4) 

{7,} {7,} <=1 

where (M) 7ii7 are the elements of the Q x Q matrix 

(e? J 1 ... l\ 



M(Q x Q) 



(5) 



V 1 1 ... 1 J 

From the structure of the matrix in Eq. (J5J), it is clear that there are only two linearly 
independent eigenvectors, so the number of nontrivial eigenvalues A and the order of the 
characteristic equation is also equal to two. This can be explicitly verified as follows. Write 
the characteristic equation |M — IA| = 0, I being the identity matrix. Successively subtract 
the second row from the first row, third row from the second, etc., until all rows have been 
accounted for. Similarly, successively add column Q with column Q — 1, Q — 1 with Q — 2 
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etc, so as to obtain a final block diagonal determinant with the characteristic equation, 

( e N _ X _ Xe N _ A 
X Q x det =0. (6) 

V 1 Q-Xj 

Neglecting the Q — 2 trivial eigenvalues, a simple change of variables A = X/Q, a = Q^ 1 
and s = (e^ J — 1)/Q yields 

detr~ A " I = A 2 - A(s + 1) + s(l - a) = 0, (7) 

V - x - A , 

which exactly coincides with the characteristic equation for the ZB model given in Eq. ($2§. 
Therefore, the Hamiltonian in Eq. (|3]) provides exactly the same thermodynamics of the ZB 
model and, hence, can be considered equivalent to it. Next, the relationship between the 
ZB and Ising models is examined. The partition function in Eq. (j4]) of the Hamiltonian in 
Eq. ([3]) can be cast in the following form, 

Q Q Q Q 

Z _ V" e^fa.iW'W.l) V" e /3J5(72, 1)5(73, ^ e /3^(7JV- 1 ,l)«5(7JV,l)_ (g) 

7i=l72=l 73=1 7iV=l 

After each of the sums, a term 5(rrii, S(ji, 1)) can be inserted, since it is equal to unity. 

m,;=0,l 

Upon changing the summation order and tracing out over 7 variables, one immediately gets 
the partition function, 

N 

Z =^ II eX P [PJ m i m i+l + ^(l - m i)] = ^ eXP [~^]' ( 9 ) 
{rrii} i=l {rrii} 

where the Hamiltonian is given by, 

- (3H zb = K zh s ^m i m i + X + // z6 ^(l -m t ). (10) 

i i 

Here the coupling K z j, = (3J represents the reduced energy of a hydrogen bond and fi z b = 
q — In (Q — 1) plays the role of a reduced chemical potential within a lattice gas formulation. 
For comparison, the Hamiltonian of a one- dimensional Ising model in a similar lattice gas 
formulation is given by, 

- f3 Hi sing = K Ising mjrrii+i + fi Ising - mj), (11) 
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where Kj sing = (3 J = K z b, and Rising — pQ- The last term on the right in Eq. (TTTT) . fii S i ng , is 
temperature dependent, whereas the analogous term on the right in f lTUj) . fi z b, is not. Herein 
constitutes a fundamental difference between the ZB and one- dimensional Ising models. 

Corroboration comes from considering the zeroes of the partition functions [4j of the 
ZB model (as defined by the Hamiltonian in Eq. (J3]), or its equivalent Hamiltonian in 
Eq. (II Op ) and comparing the zeros with those obtained for the partition function for the 
one-dimensional Ising Hamiltonian in Eq. f[TT|) . The method of zeroes of partition func- 
tions is a standard tool for identifying phase transitions in spin models 



1* 



. As the control 



variable is temperature, the zeroes must be considered in the complex temperature plane, 



known as the Fisher zeroes 



19j . In the thermodynamic limit the point where zeroes cross 



the real positive axis can be identified as the transition point. We remark, that in the ab- 
sence of long-range interactions the helix-coil transition is not a true phase transition even 
for infinite polymer lengths (N — > oo), and hence, no crossing of the real axis should be 
expected. On the other hand, recent numerical studies of Fisher zeroes for helix-coil models 
with long-range interactions, and finite chain lengths, strongly suggest a first order tran- 



sition 



20]. For short-range interactions and infinite chains, such as the case treated here, 



the distributions of Fisher zeroes are clearly different and have been discussed by Poland 
and Sheraga [4] for some limiting values of parameters. Their calculation is repeated here 
within a more general framework and results are compared with those obtained from the 
corresponding one-dimensional Ising model 19]. These different scenarios are a direct con- 
sequence of the different temperature dependence mentioned above. The Fisher zeros are 
depicted in Figure 12(a) in the Re s — Im s plane for different values of the a parameter. Here 
s = {e^ J — l)/(e 9 — 1) and o = l/(e 9 — 1) are the most convenient variables to perform this 
comparison. Poland and Sheraga in [4] have shown that at limiting values of a, Fisher zeroes 
of the ZB model lie on the unit circle, so do the edge zeroes. Expected scaling of the edge 
zeroes with a was predicted to be Re(s e d ge ) = 1 — 2a; Im(s e d g e) = 2<7 1,/2 (i — u) 1 ^ 2 , so that 
the phase transition limit is approached as a decreases, only reaching it in the o — > limit 
of infinite cooperativity. Within the microscopic formulation given by Eq. ( 1T0|) . we have 
performed a numerical check and confirmed the scaling exactly as above for the vast range 
of a values. Thus, the real part of the edge zero indicates the transition temperature and 
the imaginary part can be considered as a measure of the cooperativity or of correlations 
present in the system at the transition point 2l|. Edge zeroes of the Ising model (Fig. 12(b)), 
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as given by Eq. ( ITTj) . do not lie on the circle, clearly indicating the difference in scaling with 
a values. Both the cooperativities and the stabilities clearly differ in the two cases. While 
the transition point for the ZB model is close to s = 1, it lies in the vicinity of s = for the 
Ising model, as shown in Fig. [2]^b). As previously anticipated, a numerical simulation study 
of a three-dimensional model with long-range interactions have shown that Fisher zeroes 
occupy the same spherical region of unit radius, as in Fig. [2](a), but nearly cross the positive 
real axis unlike the short-ranged ZB counterpart j^ . 

In conclusion, we have argued that a proper microscopic description of the Zimm-Bragg 
model is not a standard one-dimensional Ising model, as often tacitly assumed in the recent 
literature, but rather a one-dimensional Potts-like model. This can indeed be shown to be 
equivalent to an Ising-like one-dimensional model having however different properties with 
respect to the usual Ising counterpart, as made evident from a comparative Fisher zeroes 
analysis. 
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FIG. 1: (COLOR ONLINE) A segment of a polypeptide chain in a trans conformation is shown. 
Parallelograms indicate the plane of a virtual peptide bond, a) Diagramatic view of a polypeptide 
chain segment where the main-chain atoms are represented as rigid peptide units, linked by virtual 
bonds through the C a atoms. Each unit has two degrees of freedom due to the rotation around 
the C a — C (torsional angle 0) and N — C a (torsional angle ip) bonds. R stands for the amino 
acid residues, all other atoms have corresponding chemical labels aside, b) Coarse-grained repre- 
sentation of a polypeptide chain: the configuration is described with the help of fictitious vectors 
li, that depend on coordinates of two neighboring C a atoms and bond angle Q{ = ir — arccos ,V +1 , 
that depends on coordinates of three carbons, and a pair of torsional angles <fii, ipi. 
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a) b) 
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FIG. 2: (COLOR ONLINE) Temperature (Fisher) zeros in terms of the s = [eP J - l)/(e q - 1) at 
different values of parameter a = 1/Q = (e q — l)^ 1 for: a) Zimm-Bragg model given by Eq. (J3j) and 
b) Ising-like model, given by Eq. (jlip . As obvious, Fisher zeroes for the a = 1/2 case are similar 
for both models, since it corresponds to the model with ^i S i ng = fi z b = field. 
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